NASA Contractor Repon 181706 


ICASE REPORT NO. 88-49 


I 


ICASE 


CONVERGENCE OF GALERKIN APPROXIMATIONS FOR OPERATOR RICCATI 
EQUATIONS — A NONLINEAR EVOLUTION EQUATION APPROACH 


I. G. Rosen 

(»ASI-C£-iei7C6) CCNViBGEKCi Ci GALERKIS tJ68-JC3e6 

AAPBCilMAlICKE FCE CEEEATCB EICCAIl 
ECUAIlCliS: A SCBiJKIAB EVCIUllCli EQUATION 

AtEECACH Final Feiort (NASA) 30 p CSCL 1EA Uaclas 

G3/64 01648C9 


Contract No. NAS 1-1 8 107 
August 1988 


INSTITUTE FOR CXHIPUTER APPLICATIONS IN SCIENCE AND ENGINEERING 
NASA Langley Research Center, Hampton, Virginia 23665 

Operated by the Universities Space Research Association 


NASA 

National Aeronautics and 
Space Administration 

Langley Research Center 

Hampton, Virginia 23665 


CONVERGENCE OF GALERKIN APPROXIMATIONS FOR 
OPERATOR RICCATI EQUATIONS - A NONLINEAR EVOLUTION 

EQUATION APPROACH 


by 


I.G. Rosen* 

Department of Mathematics 
University of Southern California 
Los Angeles, CA 90089-1113 USA 
Telephone: (213) 743-2567 


Abstract 


We develop an approximation and convergence theory for Galerkin approximations 
to infinite dimensional operator Riccati differential equations formulated in the space of 
Hilbert-Schmidt operators on a separable Hilbert space. We treat the Riccati equation as 
a nonlinear evolution equation with dynamics described by a nonlinear monotone pertur- 
bation of a strongly coercive linear operator. We prove a generic approximation result for 
quasi-autonomous nonlinear evolution systems involving accretive operators which we then 
use to demonstrate the Hilbert-Schmidt norm convergence of Galerkin approximations to 
the solution of the Riccati equation. We illustrate the application of our results in the 
context of a linear quadratic optimal control problem for a one dimensional heat equation. 
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1. Introduction 


Recently a great deal of attention has been given to the analysis of, and the develop- 
ment of approximation theory for, infinite dimensional operator Riccati differential equa- 
tions. This is due in part to the role they play in the characterization of the closed-loop 
feedback structure of the solution to linear-quadratic optimal control problems for dis- 
tributed parameter systems (see [10]). Working from Temam’s [14], [15] formulation (see 
also Barbu [3]) of a class of operator Riccati equations in the space of Hilbert-Schmidt op- 
erators on a separable Hilbert space, we develop an approximation and convergence theory 
for generic Galerkin type approximations. We make the usual assumption that our un- 
derlying Hilbert space admits a densely, continuously, and compactly embedded subspace. 
The only condition that we then require on our sequence of approximating Galerkin sub- 
spaces is the usual one that the corresponding sequence of orthogonal projections converge 
strongly to the identity in the stronger toiDology. We are able to obtain Hilbert-Schmidt 
operator norm convergence of the approximating Riccati operators, uniformly in time on 
compact time intervals. 

Our effort here differs significantly from other recent treatments of the approximation 
theory for infinite dimensional operator Riccati equations (see, for example, [6]) in that i) 
we obtain Hilbert-Schmidt as opposed to strong operator convergence, and, more impor- 
tantly, ii) our convergence theory is based directly upon, and involves only, the differential 
equation itself rather than equivalent integral equations. Of course in order to do this 
we must necessarily consider a somewhat more restrictive, but still sufficiently interesting 
from an applications point of view, class of problems. For example we require that the 
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linear part of the equation be strongly coercive on a space of Hilbert-Schmidt operators 
and that the non-homogeneous or quasi autonomous perturbation (the state penalization 
operator in the context of the LQ control problem) be Hilbert-Schmidt. Our treatment 
here is restricted to the quasi- autonomous or constant coefficient case. The temporally 
inhomogeneous problem requires a different approach; our results for the time dependent 
case will be reported elsewhere. 

Our convergence theory is based upon a generic approximation result for nonlinear 
quasi-autonomous evolution equations in Banach space with dynamics described by accre- 
tive operators. In section 2 we prove a nonliiaear analog of the well-known Trotter-Kato 
theorem on the approximation of linear semigroups - i.e. stability (uniform dissipativity) 
and consistency (strong operator convergence of the resolvent) yield convergence (see, for 
example, [11]). The result that we prove here is closely related to similar approximation 
results for nonlinear evolution systems which have appeared elsewhere in the literature (for 
example, [4] and [8]). However the theorem we prove below is one that is most appropriate 
for the particular class of problems that are of interest to us here. We follow Goldstein 
[8] and give a proof using an idea originally suggested by Kisynski [9] for the linear case 
wherein convergence is demonstrated via an application of an existence theorem to an 
evolution equation in an appropriately constructed sequence space. 

In section 3 we briefly outline Temam’s [15] formualtion of the Riccati equation as a 
well posed nonlinear evolution equation in the si^ace of Hilbert-Schmidt operators. Our 
approximation and convergence theory is developed and presented in section 4. In section 
5 we illustrate the application of our results in the context of a linear quadratic optimal 
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control problem (i.e. the linear quadratic regulator problem) for a one dimensional heat 
equation. 

2. An Abstract Approximation Result for Nonlinear Evolution Equations 

Let Xo be a real Banach space with norm denoted by | • [o- Let Xq be its dual and 

let Ao : Xo — ♦ 2^® be an, in general, multi-valued nonlinear closed accretive operator (i.e. 

{(x,y) : X £ Xo, y £ Aox) is a closed subset of XqxXo, and |xi-X 2 |o < |xi-|-Ayi-(x 2 -t-Ay 2 )|o 

for all A > 0, x; e Xo, and ji s AoXi, i = 1,2) on Xo. Define the domain of Ao to be the set 

Dom(Ao) = {x £ Xo : Aox ^ ^}, and the range of Ao to be the set 7^(Ao) = Aox. 

xeDom(Ao) 

Since Ao is accretive, it follows that for A > 0 the resolvent of Ao at A, Jo(A) = (I + AAo)“^, 
is a well defined single-valued nonexpansive operator (i.e.|Jo(A)yi — Jo(A)y 2 |o < |yi — y 2 |o, 
yi»y 2 £ Dom(Jo(A))) defined on Dom(Jo(A)) = 7^(1 -1- AAo). Suppose that T > 0, let 
t — »• fo(t) be an Xo — valued map defined on [0,T] and let Xq £ Xo. We consider the initial 
value problem in Xq given by 

(2.1) Xo(t) + AoXo(t) 3 fo(t), a.e. t£(0,T), 

(2.2) xo(0) = Xq. 

We shall say that a function Xo(-) : [0,T] Xo is a strong solution to the initial value 
problem (2.1), (2.2) if it is continuous on [0,T], Lipschitz on every compact subinterval of 
(0,T), differentiable almost everywhere on (0,T), and satisfies (2.1) and (2.2). We shall 
call Xo an integral solution of (2.1), (2.2) if it is continuous on [0,T], satisfies (2.2), and if 
the inequality 

1 1 

-|xo(t) - x|^ < -lxo(s) -^\o+ J < fo(r) - y,xo(r) - x >o dr 
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holds for all X e Dom(Ao), y e Aqx, and 0 < s < t < T where the pairing < • >o: 

Xo X Xo — » R is given by < y,x >o= sup{x*(y),x* e <^o(x)} with #o : Xq 2^® denoting 
the duality map on Xq (see [3]). A strong solution of (2.1), (2.2) is of course also an integreil 
solution. The following existence, uniqueness, and regularity results are given by Barbu 
[3]. 

Theorem 2.1 Let {X, | • |} be a real Banach space and sujjpose that A : X — > 2^ is a 
closed accretive operator on X. Let C be a closed convex cone in X such that Dom(A) C C 
and 72.(1 + AA) D C, A > 0. If x® e Dom( A) and f e Li(0, T; X) with f(t) e C, a.e. t e (0, T), 
then there exists a unique integral solution x to the initial value problem 

(2.3) x(t) + Ax(t) 9 f(t), a.e. t e (0,T), 


(2.4) x(0) = x°, 

with x(t) e Dom(A) a.e. t e (0,T). Moreover, if y is the integral solution to (2.3), (2.4) 
with f replaced by g e Li(0,T;X), then 

|x(t) - y(t)P < ]x(s) - y(s)p + 2 ^ < f(T) - g(i'), x(r) - y(r) > dr 

for 0 < s < t < T. 


Theorem 2.2 If in addition to the hypotheses and conditions of the previous theorem we 
have x*^ e Dom(A) and f e W^’^(0,T;X), then the initial value problem (2.3), (2.4) has a 
unique strong solution x with x e W^’°°(0,T;X), x(t) e Dom(A) a.e. t e (0, T) and 


|x(t)| = |Ax(t) + f(t)| < |Ax"+f(0)| +2 / |f'(s)|ds, 

Jo 


a.e. t € (0,T). 
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Furthermore, if X and X* are uniformly convex, then x is everywhere differentiable from 


the right and 

^(() = (^x(i) +/(())“, 0<(<T 

where (Ax(t) + f(t))® denotes the unique element of minimum norm in the set Ax(t) + 
f(t). 


Henceforth, we shall assume that there exists a closed convex cone Co C Xo for which 
Dom(Ao) C Co and 7?.(I + AAq) D Cq, A > 0, that x^ e Dom(Ao) and that fo e Li(0,T;Xo) 
with fo(t) e Co for a.e. t e (0, T). Consequently Theorem 2.1 applies and we are guaranteed 
that the initial value problem in Xo, (2.1), (2.2), admits a unique integral solution. 

We prove an approximation result in the spirit of the well known Trotter-Kato theorem 
for the approximation of linear semigroups (see, for example, [11]), and the approximation 
theorems given by Crandall and Pazy [4] and Goldstein [8] in a nonlinear setting. Although 
the convergence theorem we shall prove here does not differ significantly from the ones given 
in [4] and [8], the latter results are stated inappropriately and are somewhat too restrictive 
for the application we intend to consider below. Following Goldstein [8], the proof we shall 
give here is based upon an idea first suggested by Kisynski [9] in the context of linear 
semigroup approximaton. We argue convergence via an application of the existence and 
uniqueness result. Theorem 2.1. 

For each n = 1,2,... let X„ be a closed linear subspace of Xo, and let A„ : X„ — + 2^" be 
a closed accretive operator on Xn. As was the case above with Aq, for A > 0 the resolvent 
of A„, Jn(A) = (I + AA„)“^ is a well defined, single-valued, nonexpansive operator defined 
on Dom(Jn(A)) = 7^(1 -f- AAn). We assume that there exists a closed convex cone Cn in X^ 
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for which Dom(An) C Cn and 7^(1 + AAa) D Cn, A > 0. We assume further that for each 


n = l,2,...,x° is an element in Dom(An) and that f„ is a function in Li(0, T;X,i) with 
fn(t) t Cn for almost every t e (0,T). It then follows that Theorem 2.1 implies that the 
initial value problem in Xn given by 

(2.5) Xn(t) + AnXn(t) 3 fn(t), a.e. te(0,T), 

(2.6) x„(0) =xl 

admits a unique integral solution x,i e C(0,T;X„). 

In our discussions below, we shall use the notation lim Dn D Do where D„ C Xn and 

n— ►oo 

Do C Xq. By this we shall mean that for each zq e Do there exists a sequence {zn}^i with 
Zn e Dn and lim z„ = Zq. Our fundamental approximation and convergence result is given 
in the following theorem. 

Theorem 2.3 Suppose 

(i) lim Cn D Co 

n— ►oo 

(ii) ^li^fn(t) = fo(t) and there exists age Li(0, T) for which |fn(t)| < g(t), n = 1,2,..., 
for almost every t e (0, T) 

(hi) limn^ooX° = Xq 

(iv) lim Jn(A)yn = Jo(A)yo for each A > 0 whenever y„ e Cn and limn— cxsyn = Yo « Co- 

n-^oo 

Then lim Xn = Xq in C(0,T; Xq) where x„ for n = 1,2,... and Xo are the unique integral 

n — ►OO 

solutions to the initial value problems (2.5), (2.6) and (2.1), (2.2) respectively. 

Proof : If we define the linear space X over the reals by 

X = {u = {Un )^0 Un e X„,n = 0,1,2,..., lim Un = uo} 

n— ►OO 
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and for u = {un}^o ^ ^ set ju] = sup |un|o, then {X, | • |} is a real Banach space. Let C 

n 

be the closed convex cone in X given by C = {u = {un}^o ^ X ; e Cn,n = 0,1,2..} . 

Note that if Co ^ 4>, hypothesis (i) implies that C ^ <j>. Define the operator A ; X 2^ by 
' {v = {y„}^o = e .4„u„, 

Au='j n = 0,1,2, ...,lim„_oo^’n = t^o), u = {u„}^o ^ Dom(A) 

u ^ Dom(A) 

Dom(A) = {u = {un}^o e X : u„ e Dom( A„), n = 0, 1, 2 . . . , and for each n = 1,2, .. . 
there exists a v„ e AnU„ for which lim v„ = vq e AqUo}. 

11— ►OO 

The operators A„ being closed and accretive implies that the operator A is closed and 
accretive as well. Indeed, for example, for A > 0, u' e Dom(A), i = 1, 2, and v* e Au* , i = 1, 2, 


we have 

ju^ - u^l = sup \ul - U^ jo < sup ju^ + Xvl - (u^ + Av2)|o = |u^ + Av^ - (u^ + Av2)|. 

n n 

The fact that A is closed can be argued analogously. 

Clearly Dom(A) C C and it also follows that for A > 0 7?.(I + AA) e C. To see this let 
A > 0 and let v = {vn}^o ^ Then v„ e C„ and C„ C 71(1 + AA„), n = 0, 1, 2, ... imply 
that we may define Un = Jn('^)vm n = 0, 1,2, .. and set u = {un)S^o- Now v e C C X and 
hypothesis (iv) yield lim u„ = lim Jn(A)v„ = Jo(A)vo = uo; thus u e X. Also, if we set 

n— ►oo n— ^oo 

Wn = (vn-Un)/A, n = 0,1,2,.. then w = {wn}^o ^ X, %Vn e AnUn and lim Wn = wo e AoUq. 

n — ►oo 

It follows that Un e Dom(An) and therefore that u e Dom(A) with w e Au and v e (I+AA)u. 
We conclude that v e Ttil + AA) and therefore that C C 7?.(I + AA). 

Define z° = {x°}^q. Then hypothesis (iii) implies that z e X. Moreover, using (iv) 
it can be shown that Dom(A) = {u = {u„}^o e X : Un e Dom(A„),n = 0,1,2,...}. Thus 
z e Dom(A). If we define the function f by f(t) = {f„(t)}^o for almost every 1 1 (0, T), then 
hypothesis (ii) implies that f e Li(0, T;X), and fn(t) e Cn, n = 0,1,2,... a.e. t e (0, T) implies 
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that f(t) e C, a.e. t e (0, T). It follows from the auguments given above and Theorem 2.1 
that the initial value problem in X 

z + Az(t) 3 f(t), a.e. t e (0,T) 


z(0) - z® 

admits a unique integral solution z e C(0, T;X) with z(t) e Dom(A) for almost every 
t e (0, T). In addition, inspection of the proof of Theorem III. 2.1 given in [3] (paying 
particular attention to how solutions are actually constructed) reveals that we must have 
z(t) = {xn(t)}^o where for n = 1,2,..., x„ is the unique integral solution to the ini- 
tial value problem (2.5), (2.6) in X„ and Xo is the unique integral solution to the initial 
value problem (2.1), (2.2) in Xq. Consequently, since z(t) e X, t e [0, T], it follows that 
lim x„(t) = Xo(t) for each t e [0,T]. 

n — ►oo 

Finally, that the convergence of x„(t) to xo(t) is in fact uniform in t for t e [0, T] can 
be argued as it was done by Goldstein in [8]. Let e > 0 and t e [0, T] be given. Then 
z = {xn}^o ^ C(0, T;X) implies that |z(t) — z(s)| < e/3 for all 0 < s < T with |t — s| < 
26 = 2<5(e). Now lim Xn(s) = xq(s) for each s e [0,T] implies that |xn(s) — xo(s)lo < c/3 

n — Koo 

for all n > Ns = N(s, e). Thus 

Ixn(t) - Xo(t)|o < |Xn(t) - Xn(s)|o -f- |x„(s) - Xo(s)|o -|-|xo(s) - Xo(t)|o < C, 

for all 0 ^ s < T with |t — s| < 26 and n > Ng. Let k be the greatest integer less than or 
equal to T/6 and set N = max{No,N«,N 26 ,. . . ,Nka,NT|. Then |x„(t) — xo(t)|o < e for all 
> N = N(e) and the proof of the theorem is complete. 
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3. Operator Riccati Equations on Spaces of Hilbert-Schmidt Operators 

In this section we briefly review and outline Temam’s [15] results on operator Riccati 
equations set in the Hilbert space of Schmidt class operators on a separable Hilbert space. 
Let H be a real separable Hilbert space with inner product and associated induced norm 
denoted by (•,•) and [ • | respectively. Let V be another real separable Hilbert space 
with inner product < •,• > and induced norm |1 • |1 . We assume that V is densely and 
continuously embedded in H (i.e. V = H, and there exists a positive constant /x for which 
W\ ^ ^ V). Identifying H with its dual, H*, we have V <— » H = H* V* with 

the embedding of H* in V* being dense and continuous as well. Denote the usual operator 
norm on V* by (| • ||*. Let 7 denote the canonical isomorphism (Riesz map) from V onto 
V*. Then for </?, c V we have 

where (•, •) in the above expression denotes the usual extension of the H inner product to 
the duality pairing between V* and V. We assmne fm'ther that the embedding of V into 
H is compact. It then follows that 7 "^ e £(V, V*) n £(H, V) and that 7 ”^ is self-adjoint, 
positive, and compact as a mapping from H into H. We note that V* is in fact a Hilbert 
space with inner product < • >* given by 

< >*=< 7 “V. 7 ”V > = e V*. 

We have that 11<^1|* = yj< >*, v? e V*. 

Since 7 ~^ is self-adjoint, positive, and compact on H there exists an orthonormal basis, 
{ejtj^j, for H such that 7 ~^eic = k = 1,2,... for some real numbers k = 1, 2, .... 

Consequently and {pkeicj^i are orthonormal bases for V and V* respectively. 
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For separable Hilbert spaces X and Y, let HS(X,Y) denote the Hilbert space of Hilbert- 
Schmidt operators from X into Y. Denote the corresponding inner product and induced 
norm by [•,-]hs(X,Y) and | • |hs(X,y) respectively. Set H = HS(H,H) with 

OO 

= [$, ^]hs(H,H) = ^eic), 

k=l 

and l$|-w = 1$|hs(H,H) = \/[^, ^ H. Define the Hilbert space V by V = 

HS(V*,H) n HS(H, V) with inner product 

OO OO 

[$, 'Fjv = [$, 'I']hs(v,h) + [$, ^]hs(h,V) = /5k($ek, ^e^) + ^'ek >, 

k=l k=l 

ii^iiv = € V. It is not difficult to show that the dense and continuous 

inclusions HS(V*,H) C HS(H,H) C HS(V,H), and HS(H,V) C HS(H,H) C HS(H,V) 
hold. Also, it can be argued that HS(V*,H) and HS(V,H), and HS(H,V) and HS(H,V*) 
are dual pairs with respect to the duality pairing 

OO 

= ^($ek,'^ek). 
k=l 

It follows that V* = HS(V,H) + HS(H,V*) and, identifying 1-L with its dual that 
V "H V*, with the inclusions dense and continuous. 

Let a(-,-) : VxV — > Rbea bounded strongly V-elliptic bilinear form on V x 
V. That is, there exist constants a,/d > 0 for which a( 9 ,y?) > o||(^lP and |a(<^, V’)! < 
/^llv^llllV’IK 0 e V. (With minor modification our theory continues to hold if the form 
a(-,-) satisfies the coercive inequality a(<^, 9 ) + A|(,?|^ > al|<^||^, 9 c V, for some A e R. 
To keep our presentation here as simple as possible, however, in our discussions below we 
shall treat only the case A = 0.) Let ^ e £(V,V*) be the operator defined by {A<f>,‘4>) = 

10 
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a.(4>,tl>), € W. Let a*(-, •) : V x V — > R be the form which is adjoint to a(-, •); that is 

a*(^,t/>) = <i>A e V. Then a*(9,v?) > a||(^|P and \a*{<^,\p)\ < /311<^|lH'0ll, <^,0 e V. 

Let A € £(V,V*) be the operator defined by = a*(9,0) = a(0,(^), e V. 

Let Co be the closed convex cone in H given by Co = {$ e HS(H,H) :$ = $*,$> 0}. 
Let T > 0 and suppose that II§ e Co and that Fo e Li(0,T;H) with Fo(t) e Co for almost 
every t e (0,T) are given. Let $ — >• .Fo($) be a single valued map defined for each # e Co 
with range in 7i which is continuous from Ti into itself. Assume further that has the 
property that 

> 0 , 

for all $,'1' e Co. We note that if z — > fo{ 2 ) is a single valued complex function of the 
complex variable z with .Fo(O) = 0 and which is analytic on the nonnegative real axis, then 
the mapping $ satisfies .Fo(^) e for all $ e Co and is continuous from H into 

"H (see Dunford and Schwartz [5], Theorem XI.6.7.7). We seek a solution Ho to the initial 
value problem 

(3.1) rio(f) + >l*no(t) + no(t)>l + .;^o(no(t)) = Fo(t), a.e. te(0,T) 

(3.2) n„(0) = ns 

with Ilo(t) e Co, a.e. t e (0, T). 

Note that when .?^o(^) = (3.1) becomes the standard quadratic Riccati equation. 

Elementary properties of Hilbert-Schmidt operators (see [14]) can be used to argue that 
for ^ e Co we have 

= [$($ - ^»), $ - <Li]^ + [($ - ^)'I', $ - > 0. 

Essentially the same argument can be used for the case of .Fo($) = where n is any 
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positive integer. For more general quadratic terms, for example, the one that results in 
the case of the linear-quadratic optimal control problem, .Fo($) = with M t £(H), 

M = M*, and M > 0, once again a similar argument will work so long as > 0, $ e Co 
(see Proposition 2.2 in [15]). This will of course be ti'ue if $M = for all # e Co 

Remark 3.1 At this point a comment regarding the relationship between the operator 
A* and the adjoint of A is in order. Define the operator A : Dom(A) C H — ♦ H to 
be the restriction of the operator A to the set Dom(A) = e V ; A<p e H}. It can 
be shown (see [12]) that Dom(A) = H and consequently that A admits an H-adjoint, 
A* : Dom(A*) C H H. The operator A* is the extension of A* to an operator defined 
on all of V or, equivalently, A* is the restriction of A* to the H-dense subset Dom(A*) = 
€ V : A*(fi € H}. We note also that -A is the infinitesimal generator of a uniformly 
exponentially stable semigroup, {T (t) : t > 0}, of bounded linear operators on H. Similarly, 
— A* is the infinitesimal generator of the adjoint semigroup {T(t)* : t > 0} on H. In 
addition, it can be argued that both of these semigroups admit respectively restrictions 
and extensions to analytic semigroups on V and V* (see [1], [13]). 

An appropriate refomulation of (3.1), (3.2) will allow an application Theorems 2.1 and 
2.2. Define the operator Lq e £(V, V*) by Lo$ = A*^ -b for $ e V. It is not diflScult 
to argue that Lq is strongly V — elliptic; that is, there exists a constant a; > 0 for which 

(3.3) [Lo#,$]w>a>||$||2, 

and therefore that the set Dom(Lo) = {$ e V; Lq$ e W} is dense in (see [13]). Define 
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the operator Ao : Dom(Ao) CH. 'Hhy 


Ao$ = Lo^* + ^ 0(^)1 ^ e Dom(Ao) = Dom(Lo) n Cq. 

It follows that Aq is a closed operator on Ti. and that it is strongly V — monotone. That is, 
for #,’5’ e Dom(Ao) we have 

From this it can be argued at once that Aq is accretive on its domain and that 7^(I+AAo) D 
Co for all A > 0. 

We rewrite (3.1), (3.2) as the initial value problem in given by 

(3.4) rio(t) + AoIIo(t) = Fo(t), a.e. t e (0,T) 

(3.5) n„(o) = nj 

Now Dom(Ao) = Dom(Lo)nCo = Co and recall that it was assumed that Hq e Co. Also 
it was assumed that Fo e Li(0, T; 'H) with Fo(t) e Co, a.e. t e (0, T). Consequently Theorem 
2.1 yields the existence of a unique integral solution Ho to the initial value problem (3.1), 
(3.2) with Ilo(t) e Co for almost every t e (0,T). If it is further assumed that 11° e Dom(Ao) 
(i.e. that ng e Co and that X*ng + Tl^A e HS(H,H)) and that Fo e with 

Fo(t) e Co, a.e. t e (0, T), then Theorem 2.2 implies the existence of unique strong solution 
Ho e W^’°®(0,T; H) with ^*no(t) + IIo(t).4 e HS(H,H), a.e. t e (0,T). Using density it can 
be argued further that if Hg e Co and Fo e L 2 ( 0 ,T; V*) with Fo(t)* = Fo(t), Fo(t) > 0, a.e. 
t € (0, T), then there exists a unique solution Ho e L 2 ( 0 , T; V)D C(0, T;?^) nH^(0, T;V*) 
with Ilo(t) e Co for every t e [0, T] (of course in this setting, Hq is a solution of (3.4), (3.5) 
in a distributional or V* sense). 
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4. Approximation Theory 


For each n = 1,2,... let H„ be a finite dimensional subspace of H with H„ C V for all 
n. Let Pn : H — > H„ denote the orthogonal projection of H onto H„ with respect to the 
standard inner product, (•, •), on H. We shall require the assumption that 

(4.1) lim llPnVP - (/?11 = 0, e V. 

n — +00 

Note that (4.1) implies that lim |Pn<^ ~ = 0? e H and that P„ is uniformly 

n — ^oo 

bounded in n in both £(H) and £(V) with respect to the uniform operator topologies. 

It will be necessary for us to define an extension of the operator P„ to a bounded 
operator from V* into V*. We do this as follows. For ip e V* define PnC^ e H„ to be the 
representer in H„ of the functional <p restricted to a functional on Hj,. That is P„C(? = (^„ 
where is the unique element in Hn (guaranteed to exist by the Riesz representation 
theorem applied to the Hilbert space {H„, (•, •)}) which satisfies {ip,0n) = (y?n» ^n), ^ H„. 

Since Hn C V C V*, P„ can be considered to be a linear map of V* into itself. Indeed 
PnV? = is the continuous linear functional on V given by (PnV?, ^) = {‘pn, ^), ^ e V. 
The definition of the orthogonal projection of e H onto H„, (Pn P - <P,^n) = 0,0„e Hn, 
reveals that we have in fact defined an extension of P„ to an operator on V* with range 
H„ considered as a subspace of V*. Using the fact that for ip e V* cmd ip e Y we have 
(PnV?,^) = (PnP,^) - {PnP,ip - PnV’) = (P„(^, PnV’) = (<^, Pn^), and that the P„ are 
uniformly bounded in £(V), it is not difficult to show that the P„ are uniformly bounded in 
£(V*). It is worth noting that the extension of defined above agrees with the operator 
that would be obtained if we were to extend in the usual way by considering it as a 
bounded operator defined on the dense subset H of the Banach space V*. 


I 
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It follows from assumption (4.1) that lim llPnV^ ~ V^ll* = 0, e V*. Indeed, for 

n— ►©© 

9 ? e H we have ||Pn<^ ~ V^ll* ^ /^|PnV^ ~ V^l 0 as n — + oo. Since H is dense in V*, and 
the Pn are uniformly bounded in C{V*), it follows that \\Pn<f ~ S^||* — »• 0 as n — > oo for 
(f eV*. For e. V*, we have 

<Pn<^,V’>*= (Pn<^,7”V) = (V?,Pn7”V) = (V’)7”SPn7”V) = < V?, 7Pn7” V >* • 

Thus P* € £(V*), the adjoint of P„ considered as a bounded operator on V*, is given by 
P* = 7Pn7~^- Assumption (4.1) yields 

11P> - V’ll* = ll7"^P>-7"Vll = ||7 "SPi. 7~V-7”V|| = ||Pn7"V - 7"Vll 0 

as n — oo, for each <p e V*. 

For each n = 1,2,... define the finite dimensional subspace Ha of H by 

Ha = {^„P„ : e C{Ea)}. 

Note that Ha is clearly a subspace of both H and V since Hn finite dimensional implies 
that all operators in Ha are of finite rank. Define the closed convex cone Cn C Ha by 
Cn = {^nPn e Ha ‘ ^a = > 0}. Using the fact that Ha C V and C„ C Co, we 

define the operator An : Dom(An) C Ha 'Ha by 

(4.2) An($nPn) = {Ao($nPn)}|«„, ^nPn « Dom(An) = 

That is An(^>nPn) is the element Ao($nPn) in V* restricted to a linear functional on Ha- 
Since Ha is a finite dimensional Hilbert space, the Riesz representation theorem implies 
that An($nPn) = ^nPn « Ha where ^'nPn is the unique element in Ha which satisfies 

[Ao(^nPn)i ©nPn]'H ~ ['^nPiu ©iiPnj'W? ©nPn ^Ha- 
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The operators An defined in (4.2) are in effect the standard Galerkin approximation to 
Aq. Furthermore, the definition (4.2) leads to the same approximation to Aq that would 
be obtained via the more conventional procedure wherein the operator A in the definition 
of Lo is replaced by its Galerkin approximation. Indeed, for each n = 1,2,... define 
the operator An e C{Hn) by = A<^n\^ni e Hn. That is, An^n = V’n vi^here 

V>n is that element in Hn (once again whose existence and uniqueness is guaranteed by 
the Riesz representation theorem) which satisfies (>I<^n,^n) = (^n)^n)r « Hi,. Noting 

that since P„ is the orthogonal projection of H onto Hn, it is not difficult to argue that 
[#Pn,0n]-w = all $ € V* and ©„ e 'H„, it then follows that for $nPn and 

^nPn c Hn we have 

[An($nPn),^nPn]7i = [Ao($„Pn), ^nPn]>i = [Ao($uPn )Pn, ^nPn]?i 

= ^nPn “h ^nPn*^Pn *^o(^nPii)Pn? ^nPn]?i 

~ S ^nPn^kj '^nPn^k) 4" (^nPn*APnC)c , + (•^o(^nPn)Pn^k? ^nPn^k)} 

k=a 

= S {(-A*$nPnek,^nPnek) + (^^Pn^k, ^n^aPuek) + (^0 (^nPn)Pnek, ^nPnek)} 
k=l 

~ {(*Aj^$nPnGk? ^nPn^k) 4" (*^n Pn^k? ^nPn^k) 4" (Pn«^o(^nPn)Pn^k^ ^nPn^k)} 

k=l 

— 52 ({*An^n 4" ^n*An 4“ Pn»^o(^nPn)}Pn^k) ^nPn^k)j 

A;=l 

or 

An(^nPn) = + :?2i(^nPn )}Pn , 

where = Pn*^o(^)^ $ 6 Co* In particular when /o(^) = we have An{^nPn) = 

{Al^n+^nAn+^l]Pn- 

For each n = 1,2,... define II® e Cn by II® = PnlloPn and Fn e Li(0,T;?^n) by 
■Fn(0 = PnFo(i)Pn^ for almost every t e (0,T). We consider the problem of finding a 
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solution Iln to the initial value problem in Ha given by 


(4.3) ri„(t) + AnII„(t) = F„(t), a.e. t€(0,T) 

(4.4) n„(o) = n; 

with Iln(t) e Cn for almost every t c (0,T). The definition of the operator A„ together with 
the properties of the operator Ao yield that for each n = 1,2,... An is closed and strongly 
V — monotone on its domain. Hence A„ is accretive with 7?.(I + AAn) D Cn for all A > 0. 
Since Dom(An) = OT = Cn,II° e Cn, and Fn e Li(0,T;?in) with Fn(t) e Cn for almost 
every t e (0,T), Theorem 2.1 yields the existence of a unque integral solution Hn to the 
initial value problem (4.3), (4.4) with Iln(t) e Cn, a.e. t e (0, T). 

We shall argue the convergence of Hn to Ho as n — >^ oo in C(0, T; 7{) (i.e. the Hilbert 
Schmidt norm convergence of Hn(t) to Ho(t), uniformly in t for t e [0, T]) via an application 
of Theorem 2.3. In order to do this we shall require some preliminary lemmas. The first 
lemma below is a technical lemma which can also be found in [7]. For completeness we 
state it here and have included its rather brief proof. 

Lemma 4»1 If is an absolutely summable sequence of real numbers, then there 

exist sequences {bi}j^i and {cj}j^j such that lim bj = 0, {ci}j^i is absolutely summable, 

i — ►oo 

and aj = biCj, i = 1,2,... . 

OO 

Proof; Let a = ^ lajj and for j = 0,1,2,.. define the nonnegative integers kj as follows. 
i=l 

Let ko = 0 and let kj denote the first index for which 
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Set bi = 1/j and c-, = jai, for i = kj_i +l,...,kj, j = 1,2,... . Then bicj = ai, 

OO OO OO 

i = l,2,...,.lim bi = 0, and £ |cj| = 2 l^kl < « + Z) 57 < oo. 

i=i j=i k=kj_i+i j=i-' 

Lemma 4.2 Let X and Y be real separable Hilbert spaces with inner products denoted by 
< •, • >x and < •, • >Y, respectively. If # e HS(X, Y) then $ can be factored as $ = 
with e £(Y) compact and e HS(X, Y). 


Proof; Let and {yj}^x be orthonormal bases for X and Y respectively. Since 

$ e HS(X, Y) it has a representation in the form of an infinitie matrix $ [v?jj] = [< 

OO OO OO 

yi,$Xj >y] with ^ ^ (p'^j < OO. For i = 1,2,... set ai = Zv’y- Since the sequence 

t=l 3=1 j=l 

{aili^j is absolutely summable, we can apply Lemma 4.1 and obtain sequences {bi}-2i 


and {ci}i~i for which aj = bjCi,i = 1,2,..., lim bi = 0, and Z) Cj = Z ki| < 

i=i i=i 

OO 

OO. Define e £(Y) and e £(X,Y) by = Z < V^Vi >Y Vi, y ^ Y, and 

i=l 

OO OO 

= Z Z yi, X e X. Then = $, and since lim = 0 

i=l j=l OO 

/ \ 2 OO OO OO OO 

and Z ( ) = Z F Z ‘fh ~ Z = Z follows that is compact and 

i,j V'' • / i=i • ,=i • 

$2 e HS(X,Y). 


Lemma 4.3 

(a) lim lP„$Pn - = 0, ^ eU. 

n— ►OO 

(b) lim ||P„#Pn - $llv =0, $ e V. 

n — ►OO 

Proof; (a). We consider Pn to be an element in £(H). Then for # e 'H = HS(H, H) we 
have 

|P„$P„-#|« < |P„$P„-P„$|« + 1P^$-$|^ 
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< l^Pn-'^lw + |P„$-$l7i 

= |($P„r-$1w + IPn^-^lx 
= iP„$*-$*|w + 1P„$-$|W. 

Since $, $* e HS(H,H), we can apply Lemma 4.2 with X = Y = H to obtain $ = 
and = ($*)^($*)^ with e £(H) compact and c HS(H,H). Then 

\P - #1.^ = |(Pn - 1^^ 1 HS(H,//) ’ 

The fact that Pn converges strongly to the identity on H and that e £(H) is compact 
imply |Pn^ — — » 0 as n — > oo. Similar estimates for jPn^* — ^*|-H yield the desired 

result. 

(b). Let $ e V = HS(V*,H) n HS(H,V) and consider ||Pn^Pn - ^||hs(v,h)- 
Since |Pn| = 1 it follows that 

||P„$Pn - $||hS(V,H) < yPn^Pu -Pn^|lHS(V,H) + ||Pn<& - ^||hS(V,H) 

< ll^Pn - ^1 IhS(V,H) + yPn^ - ^||hS(V*,H). 

Now $ e HS(V*,H) implies that e HS(H,V*) and that ($P„)* = P^$* e HS(H,V*) 
where P* denotes the adjoint of Pn as an element of jC(V*). Then 

(4.5) \\Pn^Pn-nHS(V-,H) < ~ I1h5(H, V) + ~ $Hh 5(V ,ff)- 

An applicaiton of Lemma 4.2 with X = H and Y = V* and the strong convergence of P* 
to the identity on V* yield that the first term on the right hand side of the estimate (4.5) 
above tends to zero as n — ^ oo. Similarly, Lemma 4.2 with X = V* and Y = H and the 
strong convergence of Pn to the identity on H implies that the second term tends to zero as 
n — > oo as well. A similar argument can be used to show that ||Pn$Pn ~ ^1Ihs(H,V) 0 
as n — > oo and the lemma is proved. 
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Theorem 4.1 If Ho e C(0, T; Ti) is the unique integral solution to the initial value problem 
(3.1), (3.2) (or, equivalently, (3.4), (3.5)) and Iln e C(0, T;Wn) is the unique integral 
solution to the initial value problem (4.3), (4.4), then Iln converges to Ho in C(0, T; W) as 
n oo. That is, lim |IIn(t) — IIo(t)| « = 0 with the convergence uniform in t for t t [0,T]. 

n— ►oo 

Proof The desired result will follow immediately from Theorem 2.3 once we have verified 
that the hypotheses (i)-(iv) given in the statement of that theorem hold. If c Co then 
$n = Pn*^oPn € Cn and Lemma 4.3 implies that lim l$n“^o|'H = 0. Thus lim Cn D 

n— ►OO n--+oo 

Co. Lemma 4.3 also implies that lim |F„(t) — Fo(t)|?i = lim |PnFo(t)Pn — Folw = 0 

n— ►OO n— ►oo 

for almost every t e (0,T). Properties of Hilbert-Schmidt operators and the fact that Pn is 
an orthogonal projection yield |Fn(t) lx = lPnFo(t)Pnlx ^ lFo(t)lx c Lj(0,T), for a.e. 
t e (0, T). Consequently hypothesis (ii) is satisfied. Once again from Lemma 4.3 we obtain 
lim |n° — HqIx = lim |PnIIoPn ~ HqIx = 0* The verification of hypothesis (iv) is all 

n — ►oo n — ►oo 

that remains. 

Let A > 0 and let $n € Cn, n = 0,1,2,... with lim |^n~^o|?< = 0. Set = Jn(-^)^n 5 

n — ►oo 

n = 0,1,2,... where Jn(A) = (I + AAn)“^. Then recalling (3.3) we have 
Aa>||^n -Pn^'oPollv < A[Lo{^'„ -Pn^'oPn),’35'n -Pn^oPnlx 
— [(I "1“ AAn)^^n (1 "t" AAq)^oP n> Pn’^oP n]x 'i’oP n? P n'^^oP nJx 

+A [Lo{'Fo - Pn’J'oPn}, “ Pn^oPnlx ~ A[:F„('J'n)Pn “ .^o(^o)Pn, ’^'n “ Pn^oPnJw 

= [$n - Pn^oPn, ~ Pn'^oPnjx “ [^n “ Pn^oPn, “ Pn^oPnlx 
+A [Lo{^'o - Pn^oPnlPn, ~ Pn^oP«]x ~ A[:Fo(^„)Pn - :Fo(P„^oP„)Pn, 
^'n-PnW„]7< + A[.Fo(^'o)Pn-^o(Pn'5'oPn)Pn,^n-Pn^oPa]x 
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< [^n ~ Pn^oPn? ~ Pn’®'oPn]’« + ^ [Lo{^0 “ Pn^oPn}Pn? ’^'n “ Pn’I’oPn]-w 
+A [jFo(^o)Pn - :^o(Pn^oPn)Pn, “ Pn^'oPnlx 

< |l$„-P„$oPnllv* ||^n-Pn^oPn|lv+'^l|Lo{^'o-Pn^oPn}Pn|lv ' 

Il’i’n ~ Pn’^'oPnllv ^ l|*^o(^o)Pn ~ -^o(Pn^oPn)Pn||v» ll^n “ Pn^oPnlly • 

where in the above estimate we have used the facts that [Po(^n) ~«^o(Pn^oPn)Pn)^n ~ 
Pn^oPn] ^ 0 and Pn nonnegative self-adjoint and compact (being of finite rank) imply 
[:^o(^n)Pn - ^o(Pn’®'oPn)Pn, ’^'i. " P«^oPn]« > 0. Thus 

114'„ - P.^oPnllv < lb l*« - P»«oP„|« 

+(l/w){M*llc(V.V-) |Pnl£(H) + IMll^jv.V) |Pnl£(V) }|l^-P„^0Pnl|v 

+ 5lWo)-:^o(Pn^oPn)|« 
for some positive constant k. Now 

l$n - Pn^oPn\n ^ l^n “ ^o\n + 1^0 ~ Pn^oPn\n 

and 

lM*ll£(V,V) 1-^ «lc(H) + 1MIL(V,V) |P*nlc(V) — + l■^n|£(V)) 

which is uniformly bounded in n. The assumption that lim $n = Lemma 4.3, and 

n—^oo 

lim 1^0 — Pn^oPnl-w = 0 together with the continuity of the map 0 — > Pb(0) from H 

n— ►oo 

into itself yield lim ll^n ~ Pn’^'oPnllv = 0* Consequently the estimate 

n — +00 

lim I’i'n - ’i'olTi < K lim jj’J’n - ^ol|v 

n — ^00 n— »>oo 

< K lim 11^„ - P„^oPn|lv + K lim HPn^'oPn - ^oHv 

n— ♦oo n— ►oo 

together with Lemma 4.3 yield the desired result. 
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Remark 4.1 Although we are unable to demonstrate that the hypotheses of the existence 
result, Theorem 2.1, are satisfied when the initial value problem (3.4), (3.5) is considered 
in the space V, we can show that hypotheses (i)-(iv) of Theorem 2.3 are in fact satisfied 
in the stronger V — topology. More precisely if Hq e V and Fq e Li(0, T; V), then Lemma 
4.3 is sufficient to obtain the convergence of Iln to Ho in C(0, T; V). In the case of a linear 
dynamical system (i.e. when Ao e £(V, V*)) it is in fact possible to apply Theorem 2.1 in 
the space V and therefore obtain a complete approximation theory in V (see [1]). 

5. An Example 

We illustrate the application of our approximation results with a linear-quadratic 
optimal control problem (see [10]) involving a one dimensional heat equation (exten- 
sion to higher dimensions is straightforward). Let H = L2(0, 1) endowed with the usual 

1 

inner product, (9?, ^) = f (p(r])‘il;{r])dr). Consider the optimal control problem of finding 

0 

u e L2((0, T) X (0, 1)) which minimizes the quadratic performance index 

T 

J(u) = (^x(T,-),x(T,-)) + j (Qx(t,-),x(t,-)) +r(u(t,-),u(t,*))dt 

0 

subject to the linear distributed dynamical system 

(5-1) ^(t,??) - •^a(77)f|(t,?7) = bu(t,r^), t > 0, 0 < 77 < 1, 

(5.2) x(t,0) = 0, x(t, 1) = 0, t>0 

(5.3) x(0 ,?/) = x°(t/), 0 < < 1, 

where a e Lqo(0, 1) with &(rj) > a > 0, a.e. 77 e (0,1), b, r e R with r > 0, x° e L2(0,l), 
Q^Q ^ Co = {$ e HS(L2(0, 1), L2(0, 1)), # = > 0} Note that Q,Q e Co implies 
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that {Q<p)(ri) = / g(? 7 , and (Q<^)(»?) = f q(rj, 0<P(Od^ with g,ge L 2 (( 0 , 1) x (0, 1)), 

0 0 

9 (Vy 0 = </(6 v), (j{v^ 0 > <liv^ 0 = v)^ q{Vi 0 > 0» a.e. (77, 0 e (0, 1) x (0, 1). 

If we set V = Ho(0, 1) endowed with the usual inner product, < (p,xj; > = 

1 

J D(p(r})Di/^{Tj)dTj, and corresponding induced norm, 1| • 1|, then V* = H”^(0, 1). Define 
0 

the operator A e £(V,V*) by (A<p)(ip) = (aD(p,Di/^), v?,V’ e V. Then (A(p,(p) > 

all(7?|p,(^ e V and the restriction A of the operator A to the set Dom(A) = {c^ e V : 
Aip e H} (= H^(0, 1) n Hq( 0, 1) when a is sufficiently smooth) is positive, self-adjoint and 
the infinitesimal generator of a uniformly exponentially stable analytic (parabolic) semi- 
group of bounded, self-adjoint linear operators {T(t) : t > 0} on H. Let U = L2(0,l) 
and define the operator B e £(U, H) by {Bv){ri) = bv(»7), v e U,0 < 77 < 1, and define 

TZ e C{U) by ('lZv){r]) = tv{ri), v e U,0 < 77 < 1. 

The solution to the initial value problem (5.1) - (5.3) is given by 

t 

x(t, •) = T(t)x° + J T{t — s)Bu{s)ds, t > 0. 

0 

where for u e L2((0, T) x (0, 1)) we have used the shorthand notation t — > u(t) to denote 
the function t — »■ u(t, •) e L2(0, T;U). The solution to the optimal control problem is given 
in closed-loop, linear state feedback form (see [10]) by 

(5.4) u(t) = -n-^B*Uo{T - t)x(t, •) = -(b/r)no(T - t)x(t, •), a.e. t e (0, T), 

where Hq is the nonnegative self-adjoint solution to the Riccati differential equation initial 
value problem 

(5.5) rio(t) 4- >ino(t) -I- Iioit)A -I- ^o(t)57^■^S*^o(t) = Q a.e. t e (0,T) 
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(5.6) 


no(o) = g. 


If we set M = = (b^/r)I e £(H), then M = M*,M > 0 and, since = 

M$,$ e Co, > 0, $ e Cq. Thus Temam’s [15] theory (with Fo(t) = Q, t e [0, T], 
.T^o(^) = (b^/r)$^, and Hq = Q) presented in section 3 above yields the existence 

of a unique integral solution Ho e C(0,T; HS(H,H)) to the initial value problem (5.5), 
(5.6) with Ho e L 2 ( 0 ,T; HS(H, V) fl HS(V*, H)), Ilo e L2(0, T; HS(H, V*) + HS(V,H)) 


and Ilo(t) e Co, t e [0,Tj. (If Q is such that AQ + QA e HS(H, H) - for example 
ii g = 0 - then Ho will be a strong solution with Ho e W^’°°(0,T; HS(H, H))). Since 
7i = HS(L 2 ( 0 , 1 )),L 2 ( 0 , 1)) is isometrically isomorphic to L 2 (( 0 , 1) x (0, 1)), it follows that 
there exists a tto e C(0, T; L 2 (( 0 , 1) x (0,1))) with = 7To(t;^,7;), 7ro(t,7;,^) > 

0, a.e. (r) ^) e (0, 1) x (0, 1), t e [0, T] such that the solution to the optimal control problem, 
(5.4), is given by 

1 

u(t,)7) = -^y* 7 To(T - t; >7, ^x(t, Od^, 

0 

for almost every (t,? 7 ) e (0,T) x (0,1). 

We consider a linear spline based approximaton scheme. For each n = 2,3,... let 
Hn = where for j = l,2,...,n-l, is the jth linear spline ("hat”) function 

on [0,1] defined with respect to the uniform mesh {0, 1/n, 2/n, ..., 1}. That is 


‘fiiv) 


' 0 

1 + l 

j + 1 — ni] 

.0 


u ^ ?7 


-1 


n 


< r] < 


n 


.. < »7 < 


1 

n 


Let Pn : H —*■ H„ denote the orthogonal projection of H onto H,, with respect to the 
(•,•) inner product and define e £(Hn) to be the Galerkin approximation to A. More 


24 


precisely, we set = V’im V’n^V’n € Hn, where t/’n is the unique element in H„ which 

satisfies = (V’n^^n), « Hn. Let Qn = PnQ e £(Hn) and = P„^e£(Hn). 

Using the properties of interpolatory splines and density arguments it is not difficult to 
show that lim ||PnV? — V’H = 0, </? e Hq( 0, 1) = V and consequently that the assumption 

n — ►oo 

(4.1) is satisfied. It follows therefore, that the approximation theory developed in section 
4 yields that the solution II„ to the finite dimensional Riccati differential equation initial 
value problem 

(5.7) rin(f) + -^nnn(t) + n,i(t)^n + (b^/r)IIn(t)^ = a.e. t G (0,T), 

(5.8) n„(0) = 

with Iln(t) nonnegative and self-adjoint, satisfi.es lim !IIa(t)Pn — IIo(t)lHS(HH) = 0 

n— ►oc ’ 

uniformly in t for t e [0,T]. 

Since the basis is orthonormal, simply replacing the operators in (5.7), 

(5.8) with their matrix representations will not lead to the familiar symmetric matrix 

Riccati differential equation. For a linear operator L„ with domain and/or range in H„, 

we denote its matrix representation with respect to the basis by Ln- Define 

: [0, 1] by $n(^) = (v^i(»?), V?n”^(7))^ and set = ($„,$^) = 

f ^ai‘n)^n(i])^di]. Then = Af ^’^(aD^n, D$J), Gn = and Qn = 

0 

A4j^^(Q$n, )• I®*' = AIn^N) 2n = A4n Qn and flN(t) = AfNnN(t), then 

IIn is the solution to the initial value problem 

RN(t) H- >lJnN(t) ■+- flN(t)AN + DN(t)A4jj HN(t) = Qni f ^ (0?T) 

njv(o) = 
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with nN(t) nonnegative and self-adjoint. The approximating solution to the optimal con- 
trol problem then takes the form 

1 

Un(t,»?) = -^y 7T„(T - t;j7,^x(t,^)d^, a.e. (t,?7) € (0,T) X (0,1) 

0 

where 7r„ e C(0, T; L 2 (( 0 , 1) x (0, 1))) is given by 

for t € [0, T], (?7, e [0, 1] x [0, 1]. Our convergence theory yields lim = ttq in L2((0, 1) x 

n — ►oo 

(0, 1)) uniformly in t for t e [0,T], That is 
lim sup 

We note that some interesting open questions which are related to our effort here 
remain. For example, a Hilbert-Schmidt existence, approximation, and convergence theory 
for the case in which the operator A (and therefore the operator Ao) is time varying. The 
standard existence and approximation results in the spirit of Theorem 2. 1-2.3 for nonlinear 
systems involving time dependent operators are typically too restrictive to be of any use 
in the Hilbert space framework in which the problem would be formulated. An approach 
such as the one which was taken in [2] in the context of approximation methods for inverse 
problems for nonautonomous nonlinear distributed systems which was based upon Barbu’s 
[3] Theorem III. 4.2 may be more appropriate. 

Finally a Hilbert-Schmidt approximation theory for the steady state, or algebraic 
operator Riccati equation can also be developed (see [12]). (The algebraic Riccati equation 
arises in the context of the linear-quadratic optimal control problem on the infinite time 
interval). 
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